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INTRODUCTION 


The Space Shuttle configuration has more complex structural dynamic char- 
acteristics than previous launch vehicles primarily because of the high modal 
density at low frequencies and the high degree of coupling between the lateral 
and longitudinal motions. An accurate analytical representation of these char- 
acteristics is a primary means for treating structural dynamics problems during 
the design phase of the Shuttle program. The l/8-scale model program was 
developed to explore the adequacy of available analytical modeling technology 
and to provide the means for investigating problems which are more readily 
treated experimentally. The basic objectives of the l/8-scale model program 
are 

(1) To provide early verification of analytical modeling procedures on 
a Shuttle-like structure 

(2) To demonstrate important vehicle dynamic characteristics of a 
typical Shuttle design 

(3) To disclose any previously unanticipated structural dynamic 
characteristics 

(4) To provide for development and demonstration of cost effective 
prototype testing procedures 

This paper constitutes a progress report on the program to date. 

The work described herein has been conducted primarily under contract for 
the NASA Langley Research Center. 


DESCRIPTION OF STRUCTURAL MODEL 


The model is designed to represent the important structural dynamics char- 
acteristics of a Shuttle-like vehicle while keeping the fabrication costs low. 
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The general arrangement of the model is shown in figure 1. The original basis 
for the design was a 21.35 MN (4.8 x 106 lb) GLOW, 55.474 m (l82 ft) long 
parallel burn configuration (Grumman Design 619). Subsequently, under Rockwell 
International sponsorship, the forward solid rocket booster to external tank 
attachment design was modified to a single point connection representing the 
RIC prototype as of December 1972. 

Figure 1 illustrates a mock-up of the l/8-scale model which is approxi- 
mately 7.515 m (24 ft) from the external tank nose cone to the solid rocket 
booster tie-down plane. The total model is composed of 4 major components: 
the Orbiter, external tank, and two solid rocket boosters. The Orbiter is 
shown in figure 2, without the cargo bay doors, and in figure 3 with the cargo 
bay doors and nonstructural plastic fairings that complete the contours of the 
vehicle. 

Figure 4 shows the Orbiter fuselage under assembly and figure 5 is a 
NASTRAN plot of the finite-element model. The fuselage structural model is 
approximately 3*5^3 m (11.625 ft) long, contains 21 frame stations, and is 
constructed of 2024 aluminum. The bottom shell of the fuselage is 0.635 mm 
(0.025 in.) thick while the side walls and top shell are 0.508 mm (0.020 in.) 
thick. The cargo bay doors are made up of segments of 0.4064 mm (0.016 in.) 
aluminum sheet that are attached to the frames. The details of the attach- 
ment to the frames prevent the doors from resisting fuselage bending but 
allow them to act in resisting shear. 

The fuselage frames in the region of the cargo hay are constructed of 
aluminum sheet that has been bent to form a channel section. The tapered 
side wall channel section and the lower portion are attached back to back to 
form a U-shaped frame. The major bulkheads are of stiffened sheet 
construction. 

The delta wing shown in figure 6 consists of 6 spars and 4 ribs that are 
formed from 0.8128 mm (0.032 in.) 2024 aluminum sheet. The covers are 
0.5080 mm (0.020 in.) thick. NASTRAN plots of the finite-element model are 
shown in figures 7 and 8 . 

The fin structure, which includes only the structure from the fuselage 
to the center of gravity of the physical fin model, contains 3 spars and a 
closure rib. The webs are 0.8128 mm (0.032 in.) thick while the covers are 
O. 5 O 8 O mm (0.020 in.) thick. NASTRAN plots of the finite -element model are 
shown in figures 9 and 10 . 

A NASTRAN plot of the cargo bay doors is shown in figure 11 and a 
schematic illustrating the connection of the door shell to the door frames 
is shown in figure 12 . 
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The external tank contains four main components, the LOX tank, inter 
tank skirt, LH 2 tank, and the aft tank skirt. A NAS TRAN plot of the entire 
structure is shown in figure 13 . Tie total structure is approximately 
6.858 m (270 in.) long and has a radius of O.5O29 m (.'.9.8 in.). The orbiter 
interstage points are located totally within the LH2 part wf the external 
sank; the forward interstage at station 148.756 transmits vertical and side 
loads while the aft center-line interstage at station 21+5.7556 transmits 
thrust and side load. Inclined bars also at station 245.7556 connect with 
the orbiter at B.L. 33*75 and with the tank at B.L. 16.4631 to provide the 
necessary determinate supports. The solid rocket booster is connected to 
the tank at the forward end at the inter tank skirt. This pin connection 
transfers vertical, side, and all drag loads. The aft tank/SRB interstage 
is located at tank station 270. 988 and consists of 3 bars capable of trans- 
mitting vertical and side load as well as roll moment. 

The liquid oxygen tank, figure 14, is a shell of revolution composed of 
a conical shell, a cylindrical shell, and ttyo quasi elliptical end domes 
which are each formed from two tangential spherical segments. The overall 
length is I.98 m (78 in.). The tank is 2219 aluminum and all shell segments 
are welded at the joints. The primary gage is O.5O8 mm (0.020 in.) with the 
lower dome being 0.406 mm (0.016 in.), end the total tank structiire is con- 
nected to the inter tank skirt via a Y-ring located at the aft end of the 
cylindrical portion of the tank. 

Figure 15 shows the liquid oxygen tank connected to the inter tank 
skirt and to the forward tank/SRB interstage. Figures 16, 17, and 18 show 
additional details of the external tank structure. The LOX tank is of 
monolithic construction whereas the remainder of the external tank is of 
ring stiffened sheet construction, the sheet being thickened where large 
drag loads exist. The tank is 2024 aluminum and the overall length is 
4.27 m (168 in.). The chem-milled tank skin thickness is primarily 
0.406 mm (0.016 in.) and typically increases to 0.635 mm (0.025 in.) in 
load carrying areas such as the orbiter interstage connections. 

The solid rocket booster consists of a cylinder, a forward tank/SRB 
interstage, shown in figure 15, and an SRB aft skirt as shown in figure 19. 
The cylinder is 2024 aluminum and is approximately 3.7338 m (147 in.) long, 
5.080 mm (0.2 in.) thick, and has a radius of 0.2477 m (9*75 in.). 

A more complete description of the model design is presented in 
Reference 1. The significant structural dynamic characteristics to be 
represented in a model for various problem areas which are the basis for 
a model design are described in Reference 2. 
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ANALYTICAL MODELING PROCEDURE 



Basic Philosophy 


The entire vehicle has been analyzed by using N ASTRAL. In 
setting up the model and analysis procedures the following guide 
lines were established: 

(1) The model should be of sufficient refinement to adequately 
predict overall dynamic behavior. No attempt would be 
made to try to predict local panel motions. 

(2) The detail of modeling should be of sufficient refinement 
to allow us to predict internal load distributions that 
would be adequate for a preliminary design of the structure. 
Although we had no intention of computing internal loads we 
considered the analysis to be representative of an actual 
prototype design situation and we were interested in how 
NASTRAN would blend into a design environment. 

(5) The total structure would be analyzed by employing sub- 
structuring techniques to see how well this aspect of NASTRAN 
would blend into a design environment. NASTRAN could, in prin- 
ciple, of course handle the entire structure as a single unit, 
but we did not feel that this represented a realistic situation. 

(*0 Separate analyses of the LOX and SRB were performed to 

investigate the hydroelastic and viscoelastic capabilities 
of NASTRAN. 

(5) The NASTRAN weight analysis capability was used to calculate 
the individual component and total weights for the nonfluid 
portions of the model. A supplementary weight check was 
conducted and the NASTRAN results adjusted where necessary. 
Structural grid points were used as dynamic mass points using 
Guyan reduction as required. This procedure differs from 
Gruaman’s usual practice, vhich is to establish a weights 
model independent of the structural model. By this approach, 
unit loads on the weights model mass points are then beamed to 
appropriate structural node points. The dynamic model is then 
the same as the weights model or is a subset of it. This 
procedure inherently results in a small dynamic model and 
additional reduction schemes are not necessary. The equivalent 
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reduction takes place in the beaming of the unit loads from 
the weights model to the structural model. This method was 
not used because it would have required more extensive alters 
to the NASTRAN rigid formats, it would not use NASTRAN weight 
analysis capability, and it would have produced basic mode data 
at non physical points which might hinder test correlation. 


Overall Analysis Flow 


A schematic diagram of the analysis flow is shown in figure 20. As 
indicated the Orbiter was divided into five substructures: fuselage, 

cargo doors, fin, wing, and payload. Tne external tank was divided into 
two substructures: the LOX tank and the aft portion of the external 

tank that consisted of the inter tank skirt, LH 2 tank and aft tank 
skirt. The solid rocket booster was handled as a single unit consisting 
of the forward skirt, propellant cylinder and propellant, and the aft 
skirt. 

In the analysis each of the five Orbiter substructures was 
analyzed to produce reduced stiffness and mass matrices for selected 
dynamic points and interface attachment points. Modes for these 
components were then obtained with the interfaces held, the exception 
to this being the fuselage which was analyzed in a free-free state. 

This was done to aid in checking and to help understand the behavior 
of the combined vehicle. The five substructure stiffness and mass 
LAtrices were then merged to form total Orbiter mass and stiffness 
matrices. These matrices were again reduced by "freeing up" the sub- 
structure interface points to yield final stiffness and masu matrices 
that were used in the modal analysis. 

As mentioned earlier, seperate analyses were run on the LOX tank and the 
SKB to study the hydroelastic capability of NASTRAN and to investigate 
the effect of the viscoelastic properties of the propellant on the 
damping characteristics of the SKB. In the overall flow the SKB 
matrices * re first reduced and then merged with the Orbiter and 
external tank matrices to form a total Shuttle system of equations. 

The LOX tank was not reduced in this process. 

The aft portion of the external tank was reduced, analyzed 
seperately, and merged with the other components in forming the total 
Shuttle system of equations. 


Substructuring Procedure 


The ba 3 ic substructuring procedure for combining elements as 
presented in the NA5TRAN User's Manual has been followed with soon 
minor changes in the assumptions used, and with more extensive DMAP 
alters. These alters are written for both Rigid Format 3 which permits 
the use of more efficient eigenvalue analysis procedures 
while assembling the orbiter model and also for Rigid Format 7* The 
latter is required because the hydroelastic model of the L0~ tank 
results in nonsymmetric mass and stiffness matrices which cannot be 
treated in Rigid Format 3- The viscoelastic properties of the propellant 
also are accurately represented in Rigid Format 7. 

The analytical model is assembled in two phases. Tue flow diagram 
for the analysis is shown in figure 21. In the first phase, each sub- 
structure is analyzed and checked separately. The output from this 
phase is assembled onto a copy tape for the symmetric and anti- 
syranetric cas-'s and then coupled in Phase 2. 

The following changes to the basic substructuring assumptions have 
been made in formulating this procedure: 

Any external supports present are included in the Analysis 
Set (a-set). 

Any zero stiffness degrees of freedom and symmetric or anti- 
symmetric boundary constraints at the model plane of symmetry, 
are included in the 3ingle Point Constraint Set.. No other 
degrees of freedom are included in this set. 

Massei which are associated with zero stiffness degrees of 
freedom will be lost unless these degrees of freedom are 
"beamed" to adjacent points using Multipoint Constraints. 

The interface degrees of freedom may be sequenced differently 
and in different coordinate systems in any two substructures 
to be coupled. Multipoint Constraints are used to relate the 
appropriate degrees of frew-dom irrespective of local coordinate 
systems or initial sequencing. 

although the general theory presented in the NASTRAN User's Manual for 
substructuring is correct, it does not provide analysis checks at various 
critical points in the procedure. Structural plots provide analysis 
checks in this substructuring procedure but are not considered sufficient 
for verifying more than structural topology. 


The following checks have been incorporated in the analysis by 
means of extensive DMAP alters: 

A rigid body check is made in Phase 1 afte_- the generation 
of the reduced stiffness and mass matrices. Temporary 
rigid body suppo* cs are included in the deck as SUPORT 
cards for this purpose. 

The structural transformation matrices Gm, Go and D are 
used to generate equilibrium matrices for the various 
constraint ~ets ercept single point constraints. These 
equilibrium matrices represent resultants about a chosen 
origin due to unit applied loads at the appropriate 
degrees of freedom. 

Provision is- made to compute either free-free modes or 
free modes with the substructure held at th - ' interface. 

ThiB is necessary if each substructure is to be checked 
independently in Phase 1. 

A rigid body mass matrix relative to the basic origin is 
computed and compared with the general mass matrix 
calculated by the Grid Point Weight Generator. This check 
verifies that no mass lias been lost in the reduction process. 

The DMAP statements to perform these functions for Rigid Format 3 
are presented in the Appendix. 


Finite Element Model 


The number of grid points and elements used in the five (K iter 
substructures are shown in table 1. The fuselage shell structure v zz 
modeled uoing CQDMEM2 elements, a new element in NASTRAN but one that 
has been used widely at G. unman. It is essentially a quadrilateral that 
is composed of four triangles which have a common central node defined 
by the intersection of lines that connect the midpoints of the opposite 
sides of the quadrilateral. The four corner nodes need not lie in a 
plane. The fuselage U frames (see figures 4 and 3) and keel were 
Idealised using CRQD and C SHEAR elements. Here effective cap areas 
were calculated for the CRQD elements to represent the appropriate 
bending behavior. CBAR elements with appropriate offsets were used to 
represent thin ring type frames such as the engine compartment closure 
frame (see figure 2). 
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The webs in the wing ribs and spars (see figures 6, J, and 8) 
were idealized with CSHEAR elements. Again the effective rib and 
spar bending material was incorporated into CROD elements 
in the upper and lower covers. The covers themselves were represented 
by CQDMEM2 elements and some CTRMEM elements that occur at the leading 
edge. Intermediate node lines that lie between the spar and rib node 
lines were established to further refine the grid used in the covers. 

The geometry of these lines is essentially set by the location of 
fuselage frames. 

The idealization of the fin follows closely the same scheme used 
in the wing (see figures 9 and 10 ). 

The shell portion of the cargo bay door (figure ll) was idealized 
using CQDMEM2 elements with the exception of a few CQUAD2 elements 
that were required for local stability to provide an attachment point 
of the doors to the fuselage. Tne door frames were idealized as 
CSHEAR and CROD elements. Note that these frames contain two webs 
(figure 12), one common lower cap, and two upper caps that connected 
(( to the ( forward and aft shell segments. This allows the doors to 
breath in longitudinal direction. 

Although provision was made for testing four payload configurations, 
the analysis included only one that represented the full up payload of 
289 kN (63 000 lb). The stiffened box section payload was represented 
by a series of CBAR elements. The payload is shown mounted in the fuse- 
lage in figure 2. 

The fluid in the LOX tank was represented by a network of four 
concentric fluid rings, 13 levels deep. The shell was idealized as 
CQUAD2 and CTRIA2 plates while the Y-ring was represented by CBAR 
elements. The shell was divided into 22?° segments in the circumferential 
direction and IT stations in the meridional direction. 

The aft portion of the external tank (figure 13 ) was modeled using 
CQUAD2 elements to represent the shell. Double frames exist at the 
f or> ard and aft portions of the inter tank skirt and an additional 
longitudinal nods line is picked up in this region to account for the 
SRB drag attachment and the stiffening that exists in the shell. Five 
heavy frames exist in the eft external tank; the first at STA 99.98 
which is The forward tank/SHB interstage; the second at STA 148.756 which 
is the orbiter forward interstage; the third and fourth at stations 
229.156 and 245.7536 which pick up the orbiter aft interstage fitting; 
and the fifth at station 270 . 988 which is the aft tank/SRB interstage. 
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These heavy frames have internal struts to provide additional stiffening 
to the interstage attachment points (figure 18). The remainder of the 
frames are light and are included to prevent shell buckling. In the 
prototype design, the real shell was stiffened in the longitudinal 
direction. In the model this stiffening plus the skin thickness was 
lumped to yield an effective thickness which was then scaled to the 
dimensions of the l/8-scale model. This was done for the sake of 
economy in constructing the l/8-scale model. 

The solid rocket booster finite element idealization consists of 
CQUAD2 plate elements (containing membrane and bending properties) to 
represent the skin, straps, and plates; three-dimensional elements to 
represent the propellant; and offset bar elements to represent tne frames 
and longerons. A NASTRAN generated plot of th( outer shell is shown in 
figure 22 along with the frame stations. The thickness of the forward 
skirt varies from 1 to 6 mm (0.040 to 0.230 in.), the propellant cylinder 
thickness is 5 mm ( 0*1875 in.) and the aft skirt thickness is 2 mm 
(0.062 in.). The propellant is modeled by three layers (in the radial 
direction) of three-dimensional elements whose properties are 

E 1 = 172.37 MN/m 2 (25 x 103 psi), v = Oj+9, p = 1716.15 kg/m3 
(0.062 lb/in3) and a structural damping factor p = O.52 where 
3 = G'Vg' = E"/E' (E = E' + Ell, G = G' + G"). The total weight of 
the structure and propellant is 11 kN (2520 lb) . 


ANALYTICAL RESULTS 


Orbiter Component Analysis 


The analysis of the separate components conducted as part of 
Phase 1 is used to establish confidence ir the finite element models 
at that level. The NASTRAN generated weights were compared with those 
determined independently and discrepancies were rectified. The vibration 
eigenvalues and eigenvectors were calculated for the components restrained 
at their supports, or free, whichever seemed most applicable. These vere 
examined and any departure from anticipated results was investigated. 

This check helped uncover problems in the way constraints were specified 
and some other data difficulties. The lowest frequency modes obtained 
during these component analyses were as follows: 
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Fuselage : 


Free, symmetric 

Free, antisymmetric 

. . . . 62.2 Hz 

. . . . 89.I Hz 

129.9 Hz 
128.3 Hz 


Payload : 

Restrained, symmetric 

Restrained, antisymmetric .... 

. . . . 81.2 Hz 

. . . . 68.6 Hz 

268.5 Hz 
175.4 Hz 

627.7 Hz 

462.8 Hz 

Cargo doors: 

Free, symmetric 

Restrained, antisymmetric .... 

. . . . 4.6 Hz 

. . . . 156.4 Hz 

10.7 Hz 
622.2 Hz 

17.6 Hz 
1054.6 Hz 

Wing: 

Restrained 

• • • * 7T • ^ Hz 

158.3 Hz 

259.9 Hz 

Fin: 

Restrained, symmetric 

Restrained, antisymmetric .... 

. . . . 264.2 Hz 
. . . . 107.8 Hz 

841.3 Hz 
407.2 Hz 

1263.3 Hz 
1018.7 Hz 


Total Orbiter Analysis 


After the individual components were analyzed, the entire orbiter 
vehicle was coupled and a vibration analysis was performed in Rigid Format 
3. PIOTEL elements were used to connect the grid points retained for 
plotting purposes. In order to examine the deformation more readily, 
both a side view and a bottom view were plotted for each mode. Only the 
latter includes the payload. The deformed shape wa3 plotted together 
with the X, Y, and Z vectors from the underformed location. The first 
two elastic modes are shown in figures 23 to 26. The first mode at 
53*0 Hz exhibits fuselage vertical bending, fin pitching, and wing 
motion. Wing motion appears to be due to flexibility in the root 
restraint and the deformed shape is almost a straight line. The 
maximum motion point is at the fin tip and results from pitching of 
the back part of the model. The second elastic mode at 62.6 Hz is 
principally wing bending with some payload and fuselage vertical bending. 

Initial comparisons with test data indicate that there is more 
flexibility in the fin and wing attachment in the physical model than 
was allowed for in the analyses. The orbiter finite element model is 
readily adapted to exploring these effects and several runs were made 
varying the fin attachment. Results showed that the aft frame in the 
orbiter offers little stiffness to the aft fin spar in the vertical 
direction, but the forward frames are very significant. The first 
symmetric mode calculated with the forward frame vertica „ forces 
eliminated from the NASTRAH model is shown in figures 27 and 28. The 
frequency dropped from 53.0 to 48.0 and the relative deformation of the 
fin iB easily noted. 
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LO^ Tank Analysis 


The full l £>2 tank model, with no omits, was analyzed for the 
zeroth and first harmonics only for the frequency range from 8 riz to 
135 Hz. This frequency range was selected to avoid calculating the 
slosh modes which were not considered significant in our analysis. 

The modes obtained can be most readily characterized by the variation 
in pressure. 


SUMMARY OF HYDROELASTIC MODES 


fl/8-scale LOc tank! 


Frequency, Hz 

Characteristic pressure pattern 

Zeroth pressure harmonic (circumferential pressure = cos 00) 

22.9 

No nodal surfaces 

75.2 

1 node at about midtank 

91.5 

2 nodes 

115-2 

3 nodes 

First pressure harmonic (circumferential pressure = cos 19) 1 

19.2 

No nodal surfaces 

60.5 

1 node 

110 

2 nodes 

134.3 

3 nodes 


The corresponding grid point deformation for the original 
structural idealization indicated irregularities associated with the 
fin5-te element model of the lower dome. Since the pressure gradations 
in the lower hydroelastic modes were relatively uniform, it appears 
suitable to investigate the effects of dome finite element size and 
geometry using static pressure loading to save computer time. The 
static loading produced deformations very similar to those In the 
fundamental hydroelastic modes. One modification attempted, the use 
of membrane elements in place of plate elements, gave no appreciable 
improvement. The original finite element grid was then refined by 
adding more elements, and the geometry was corrected. The resulting 
deformation pattern was considered acceptable. The current version of 
the tank dome finite element represent.. t ion is shown in figure 29. 

Both the undeformed shape and the pattern under a uniform pressure of 
6.9 kN/m? (1 psi) are shown. 
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External Tank Analysis 


In assembling the external tank model, the LH 2 tank including the skirts 
at both ends was analyzed as an empty free shell. Vibration modes resulting 
from these computations indicated that above the first bending mode at 
139.2 Hz, the modes of the central portion of the LH 2 tank in the areas of the 
light frames exhibited radial deformation typical of shell modes in cylinders. 

An interesting comparison of the NASTRAN calculated weights and those 
determined independently by a weights engineer is as follows: 


Tank weight as calc -dated from structural drawings, 

including fittings, fasteners, etc 603-2 N ( 135 - S lb) 

Finite-element model weight (twice the half tank 

weight) 589.4 N (132.5 lb) 

c.g. position aft of forward dome as calculated by 

weights engineer 1.901 m (75-14 in.) 

c.g. position as calculated by NASTRAN for the finite 
element model 1.912 m (75*29 in.) 


The weight of the LH2 is distributed as nonstructural mass in the CQ,UAD2 and 
CTRTA2 elements. 

After the LH2 tank model is checked, it is reduced and coupled with the 
I&2 model. Analysis for this coupled structure has not yet been completed. 


SRB Analysis 


In order to obtain a guide for the finite element idealization of an empty 
tank, the SRB was modeled as a cylinder of radius 0.25 m (10 in.) and length 
5.O8 ai (200 in.). The finite element idealization consisted of 21 bays along 
the length and 12 bays around the circumference. The following table repre- 
sents a comparison of results between NASTRAN using the Givens method, 

Grunman’s STARS-2V program, and NASA Langley's SRA program (refs. 3 and 4, 
respectively) . The STARS-2V and S 1 programs are based on thin-shell ortho- 
tropic theory. The accuracy of the NASTRAN results are relatively good for the 
lover modes and depend upon the relative complexity of the eigenvectors. 
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EMPTY CYLINDER VIBRATION ANALYSIS 


Frequency, Hz 


— 

STARS -2V 

SRA 

NASTRAN (householder 
method) 

$ Error 

52.0 (n=2, 1st) 

51.56 (n=2, 1st) 

55.2 

6 

52.4 (n-2, 2nd) 

51.66 (n=2, 2nd) 

54.9 

5 

66.6 (n=2, 3rd) 

66.04 (n=2, 3rd) 

73.9 

11 

119.3 (n=l, 1st) 

120.-':6 (n=l, 1st) 

122.5 

3 

120.4 (n=2, 4th) 

— 

171.8 

42 

147.1 (n=3, 1st) 

— 

165.1 

12 


The undamped vibrational inodes for the full cylinders are tabulated in the 
tables that follow. The modes of most interest are the 1st and 2nd bending 
modes and the longitudinal rod and thickness shear mode. Figures 30(a) and 
30(b) show cross sectional views of the vibrational motion, and figures 31(a) , 
31(b), and 31(c) show orthographic views of the motion obtained from the 
NAS TRAN analysis. The first table also includes the results for simple beam 
theory for the modes of interest (bending and longitudinal ) based ou the com- 
posite properties 0 ? the tank. Using a structural damping factor of O.52 for 
the propellant elements, the complex eigenvalues for the lowest bending and 
longitudinal modes were obtained (Rigid Format 7) and compared with the 
undamped modes as tabulated in the second table. Simple beam theory (no shear) 
predicts a value of l/Q = 0.028, which agrees with the bending mode. The dif- 
ference between this value and that for the longitudinal mode is due to the 
thickness shear effects. (See figure 30(b).) It was found that the damped 
vibrational analysis was run more efficiently by analysing the undamped system 
first in order to narrow the search range. However, computer ruining times 
were still quite long. 
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VIBRATION ANALYSIS OF FULL PROPELLANT CYLINDER 


Mode 

Frequency, Hz 

NASTRAN 

Simple beam theory 

n=l, m=l 

56.4 

58.4 

n=0, torsion 

171.4 


n=l, m=2 

173.0 

161.0 

n=0, longitudinal 

196.1 

100.2 


VIBRATION ANALYSIS USING DAMPED SOLID FINITE ELEMENTS 



Frequency, Hz 

Damping value, 

l/Q 

(a) 

Mode 

Undamped 

Damped 

Bending - 1st 

56.38 

56.39 

0.027 

Longitudinal - 1st j 

196.0 

197.1 

.056 


a l/Q = J\ where T\ is the equivalent damping constant; 
c.f., Tong, Kin N. : Theory of Mechanical Vibrations. John 

Wiley & Sons, Inc., I960, p. 15. 


Total Vehicle Analysis 

At the time of this writing, vibration analysis results for the 
completely coupled shuttle configuration were not available. 




NASTRAN EXPERIENCES 


Hydroelastic Analysis 


Some difficulties were encountered in attempting to run the 
hydroelastic analysis. Eefore setting up the l/8-scale model LO 2 
tank, the program was run for a small problem containing 86 degrees 
of freedom in the analysis a-set. After this had been run 
successfully, the LOo tank which had 7^7 a-set degrees of freedom, 
was modeled and submitted for computation. A summary of the 
difficulties encountered . conducting the larger hydroelastic 
analysis are as follows: 

(a) Hydroelastic problems will not run in Level 15-5 of NASTRAN. 

A system 0C1 error occur 3 while executing module GKAD- This 
error has been reported to NSMO and is listed as SPR . 

(b) Often only a single Eigenvalue is extracted, using the Inverse 
Power Method, although more are present. This we now feel is 
a function of incorrect completion codes. This error is now 
listed as SPR 995 . 

(c) Fluid rings must be input in ascending order on RINGFL cards 

or program terminates with error No. 2001. This error has been 

reported to NSMO and is listed as SPR 1017* 

\ 

(d) Fluid element identification numbers are limited in size to 
approximately 30 000 or less. Numbers greater than this 
cause a 0C5 system error in Module TAl. This error is now 
listed as SPR 1016. 

(e) BAR0R card causes fatal error in hydroelastic analysis. This 
error has been reported to NSMO- 

(f) Data block Maa is not pooled correctly in module SMP2 in 
Level 15 .I. This causes fatal error 1105 if program is 
checkpointed. (Problem runs without checkpoint). This 
error has been reported to NSMO. 


One continued difficulty was the large amount of computer r mining time 
required for the eigenvalue solutions in Rigid Format 7* 
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No information is available either in the literature or from NSMO 
regarding the reduction of the number of D.O.F. when using fluid elements in 
hydroelastic problem. And yet, if shuttle hydroelastic analyses are 
to be accomplished in moderate computer time then a major reduction 
seems advisable. In order to determine i: ha reduction is 

possible, a small hydroelastic problem vat- used. It was found that 
using the internally generated fluid p int numbers on OMIT cards 
did not violate any NASTRAN rules and the program ran successfully to 
completion. These internal numbers may be calculated following the 
rules in the NASTRAN User* s Manual or an unreduced problem may be run as 
far as GP4 with diagnostic 21 turned on. 

A review of the frequencies shown in the following table indicates that 
the results are comparable for the lower frequency modes. 


EFFECT OF REDUCING FLUID POINTS IN HYDROELASTIC ANALYSIS 

[simple l/8-segment of hemispherical tank, 
total degrees of freedom = 150 



Frequency, hz 

Mode 

No omitted points 
analysis D.O.F. = 85 

Omitted fluid points 
analysis D.O.F. = 77 

1 

283 

292 

2 

421 

436 

5 

536 

544 

4 

606 

698 

5 

6 97 

797 
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SRB Analysis 


For the empty propellant cylinder the inverse power method found errone- 
ous roots and left out some roots. These roots were subsenuently found using 
the Givens method and the erroneous roots did not appear. The Givens method 
generally did not work for large problems on level 15*5 but did work on level 
15.1. The damped vibrational analysis, using Rigid Format 7, gave a fatal 
error message after finding the eigenvalues. The eigenvalue running times 
tended to be long (1000 CPU seconds on an IBM 370/165 for 176 reduced D.O.F. 
using the DFT method). These errors did not occur for very small prototype 
problems. Other difficulties that were encountered included erroneous fatal 
messages} for example, a U602 message was encountered for a singular matrix. 
These errors also seemed to be a function of the large size problems under 
consideration. 


Orbiter Coupling Analysis 


Once the EMAP alters were debugged, essentially no major problems 
were encountered as far as obtaining results for the orbiter. The 
incorporated checks and plots proved to be major aides in "debugging" 
the input data to Phase 1. Experience with the various alters is 
listed below: 

(1) Incorporating the rigid body checks in phase 1 is essential 
in determining if there are any erroneous constraints in the 
substructures . 

(2) If the rigid body check is not satisfactory and the erroneous 
constraint is limitted to a single constraint, then printing 
the reduced rigid body support stiffness [xj and obtaining 
from it the resultants of the rigid body forces helps in locating 
the coordinates of the erroneous constraint. 

(3) If the trouble is caused ty MPC's then the resultants of the 
m-set loads helps in locating MPC errors. 

(4) If MPC's and SPC's are in error, mode plots are helpful in 
locating erroneous CPC's and sometimes MPC's. 
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(5) If free-free modes are obtained, then printing the member 
forces and/or SPC forces for the rigid body modes irn.y 
help in locating the erroneous constraint since the 
structure should be free of stress. 

(6) Mode plots in phase 1 have helped in determining whether the 
appropriate nodes have been selected as dynamic degrees of 
freedom. In some cases "soft spots" were accidentally 
selected for the a-set and these caused local motions to 
show up in the mode plots. 


In order to obtain plots in the coupling run (phase 2) it is 
suggested that grid points rather than scalar points be used in the 
coupling phase. PLOTEL elements were then used to connect; the grid 
points creating a pseudostructure that is suitable for plotting. The 
grid points established in phase 2 were the grid points that are 
associated with the substructure a-set degrees of freedom. All non- 
strainable D.O.F. were removed by SPC's. It should be noted that 
each substructure had a unique grid point numbering system so tint 
the grid cards in the a-set of each substructure could he duplicated 
and incorporated in phase 2. Common interface points were made common 
by MPC's. 

If necessary the a-set of a given substructure was increased so 
that a more realistic plot could be obtained. This also necessitated 
having x, y, and z D.O.F. at all points to .->e plotted so that all 
significant motion is displayed. 

To prevent loss of mass, it is recommended that mass should not 
be assigned to grid points having nonstrainable D.O.F., such as, 
intermediate grid points in a planar frame. If assigning mass to 
such nodes is necessary, then MPC's should be used instead of SPC's 
to remove the singularity from the stiffness matrix; this will conserve 
the total mass distribution. 
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System Experience 


One of the fallouts of our analysis of the l/8-scale mod'll has 
been a further evaluation and demonstration of the program that is 
scheduled to eventually replace our own in-house system. Partly as 
a result of this work, we believe that NASTRAN is ready to handle the 
analysis of large aerospace vehicles such as the shuttle. We would 
like to point out, however, some additional features associated with 
NASTRAN which must be given consideration. 

(1) The learning curve for NASTRAN is rather flat. If you 
want to be in a position of making extensive alters to 
the rigid formats, and any aerospace company faced with 
large complex problems must be in this position, then the 
investment in learning time is large. Future levels of 
NASTRAN should concentrate on building a system that is 
more easily altered. We feel that it is far more important 
to devote NASTRAN funds to developing a sound basic system 
than to adding capability for solving specialized problems. 

(2) Our in-house developed postprocessor for converting selected 
NASTRAN element corner forces, for example, membrane elements 
and rod-shear panel assemblies, has been completed (available 
from level I5.5). Although not used on the l/8-scale model 
analysis this program is a necessity if we are to obtain internal 
member loads in a form required by our designers. 

(3) Experience in running large problems in NASTRAN should be 
established prior to actual run submissions. Adequate time 
should be provided for difficulties encountered the first few 
times a new problem i . run. The availability of experienced 
computer systems anily. ~s capable of assisting in such 
difficulties helps ma uiially n. expediting NASTRAN analyses. 
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APPENDIX 


NASTRAN SUBSTRUCTURING ANALYSIS FOR NORMAL MODES 
ALTERED RIGID FORMAT 3 FOR PHASE 1 OR 2 


Incorporated New Bulk Parameters 


(1) NOSUB 

(2) TPCOPiao 

(3) TPNAME 

(4) RMODESO 


(5) TPNAME9 

(6) TPCOmaC 

( 7 ) TPNAMEN 


Number of substructures to be coupled in this run. 
Default * 1, which indicates a phase 1 run, where 
one substructure will be reduced. 

Will put reduced stiffness and mass matrix (Kaa & 

Maa) on tape PTT. Defcul* = -1 

Label name of INPE. Use only when TPCOPfaO 

Causes restrained free modes to be obtained. The 
te&iraints are defined in an input column partition 
matrix (CPAJC), which will partition the a-set into 
J 4 C sets. Default = -1. In this case free-free 
modes will be obtained if there i3 a SUPORT card in 
the BU1K data, defining the rigid body supports. 

Although (CPAJC) is not used when FMODE * -1, it 
mist be defined in the BU1K data. It is sufficient 
to defint it as a 1 y 1 matrix. Also, don't forget 
the EIGR card if modes are to be obtained. 

Label name of INP9, which contains the column 
partition vector, reduced stiffness, and mass for 
each reduced substructure. The column partition 
vectors are used to merge the reduced 3tiffness and 
mass of each reduued substructure into a common pseudo- 
structure lineup. Use this parameter only when NOSUB>« . 

Will put th«* pseudostructure eigenvalues and eigen- 
vectors in substructure lineups on tape (INP1, INP2, 
etc.) for further processing, in this case final 
substructure mode shapes. Default ■ -1. 

Common label name of IHPl, INP2, etc. Use cnly w'. -„n 
TPCOPYNaO. 
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Incorporated New Bulk Input Matrices 


The following input matrices must be defined in the BUIK data 
on DMI cards whether they are needed or not. If they are not needed, 
defining them as a 1 x 1 matrix will suffice. 

(1) EQR Matrix 

This matrix expresses the resultants about an origin, 
due to unit rigid body support loads. The rigid body degrees 
of freedom are defined on the standard NASTRAN SUFOET card. 
The EQR matrix is necessary if the checks, which are 
incorporated in the ALTERS, are to be performed. The origin 
chosen, should be the same origin defined on the standard 
GRDPNT parameter card. 

(2) CPAJC Matrix 

This matrix is used when restrained-free modes are to be 
obtained (KMODE=l) . This matrix is a column partitioning 
vector which defines the restrained degrees of freedom from 
the analysis set (a-set) degrees of freedom. 


IMPORTANT NOTE : 

When doing a coupling run, where all substructures have been 
reduced and on tape, it was necessary to input in the BUUC data at 
least one element, to prevent a fatal error in module TA1. A thin, 
string-like rod will suffice. The element must be counted as a 
substructure so that the new NOSUB parameter was increased by one. 


Alters Incorporated (General Flow) 


STAETI 


.... f ... 


ALTER 2 

... 


r 


next 

page) 


Define 

Define 


new parameter defaults 
parameter TRUE = -1 


r 

] 

i 

-l 


NOSUB = 1 
TPCOPY = -1 
TPCOFIN = -1 
RMDDE = 1 


a 


( 


V 


w 
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ALTER 23, 23 1 

Page 193 


ALTER 24 |— 

=£=; 
ALTER 23 \ ~ 


ALTER 27~}- 
ALTER 29 |- 


ALTER 30 | 

— 

ALTER 37, 37 
Page 195 


l 

ALTER 40,4o}- 


ALTER 48, 481 



ALTER 49 f 

\ 

> 

ALE 

s 

CO 

>* 

CO 

1 

1 

1 

1 

1 

1 

1 

1 


Delete original error (possible not to have 
structural elements in BUIK for coupling run). 
Mso defines coupling parameter and other 
parameters and performs appropriate purges. 


Insert LABEL 124 (referenced in ALTER 23, 23) 


Insert COND 127, NOSIMP (Skips generation of 
stiffness Kggx if no structural elements) 


Insert LABEL 127 (referenced in ALTER 25) 


Insert 


COND 129, COUPLE 
' JUMP L30 
LABEL L29 


Possible not to have mass 
in BUIiC for coupling run. 


Insert LABEL L3G (referenced in ALTER 30) 


If this is a coupling run (phase 2), this alter 
will read in reduced substructure matrices from 
tape and merge them into pseudostructure g-lineup. 
The combined stiffness and mass are then added to 
the unreduced substructure stiffness and mass 
(if any). 

Change input Kggx to Kggy in SIjA 3 module (Kggy 
includes reduced substructure stiffnesses, if any) 


Change Mgg to Mggy in EQUTV statement. Mggy 
includes reduced substructure masses, if any. 
Also changes the mass matrix Mgg to a weight 
matrix and prints it out (wgg = Mgg x 386.4) 


Insert 


COND L49. COUPLE 
“ JUMP LBL4 
, LABEL L49 


Skips GPSP module 
if a coupling run 


Change Mgg to Mggy in MCE2 Module (Mggy includes 
any reduced substructure mass if any.) 










SEEMAI KAA, MAA (pictoral view to see if there 
are zeros on diagonal 



This package can be considered a quality control 
module. It generates checks to ensure that 
there are no errant constraints and no mass 
has been lost. See detailed flow. This alter 
also puts the reduced stiffness and mass on tape 
provided the parameter TPCOPY is not less than 
zero. 

Obtains restrained free modes (fixed at interface) 
or directs flow to obtain free-free modes. 


Change MI to M4 in original statement 


Insert LABEL L93 (referenced in ALTER 89 , 89 ) 


If a coupling run (phase 3) , this alter extracts 
out eigenvalues and eigenvectors for each reduced 
substructure and puts them on separate tapes. 
These tapes can then be used in a PHASE 3 run 
to obtain detailed substructure mode shapes. 
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•ff jNp 1 *.-. **,- 


ALTER 23 , C 3 DETATTJID FLOW 


ALTER 


\ Y 

DEFINE COUPLING PARAMETER 

) * 

COUPLE = NOSUB-2 


- If -I, this is not a 
coupling run (Phase l) 
This run would strictly 
be a Phase 1 run, when 
one substructure will be 
reduced. 


^ § 


NOEIKT = 

YES / , X NO r 

NOEIKT = +1 

f \nosihp yS 



'No elements 
defined in BUIK. 

> - m 




- v 


: v 
% 



No structural elements in BUIK. 


LABEL ER23 

PR INT PARAMETER 
NOINP 


1 //'--NO . ■■■<« 1 

t /is Nr ^no 

PURGE & CHKFNT 

OGPST 



Poaaible Cause of Error 
If coupling run, the NOSUB 
parameter must be Input in 
BUIK. 



( E ND OF ALTER 23,23*) 
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AUER 37, 37 DETAIIED FLOW 



JUMP LOOP3T c~ Top of Loop which combines 
LABEL LOOP3T reduced Substructure matrices 
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MERGE SUBSTRUCTURE REDUCED STIFFNESS & MASS 

INTO PSEUDOSTRUCTURE G-LINEUP, USING COLUMN 

PARTITIONING VECTOR CPGi 


^TmtI 


NOT 1st PASS 


LABEL L37A 
^^ASSl^ 

X<o/ 


YES 1st PASS 


ADD ABOVE MATRICES TO PREVIOUS ACCUMULATED 
SUBSTRUCTURE MATRICES 

[Kggs] < B ^ UIV - [Kgt'i = [Kggs] + [Kggil 
EQUIV 

[Mggs] [Met} = [Megs'] + [Mggi] 


LABEL L37B 


DEFINE SKIP LOOP PARAMETER 
SKIP2 = NOSUB-PASS 


INCREASE LOOP PARAMETER 
PASS = PASS + I 


/ IS \ 

'noelmt 

Vo/ 


Elements are present 
in BULK 


No elements present 
in BUIK. All substructures 
have been reduced and on 
tape. 


LABEL L37C 

yisV 

/SKIP 2> 

\<o/ 


No more reduced 

substructures on tape.- 


REDEFINE SKIP LOOP PARAMETER 
SKIP2 = SKIP2 * 1 


Although this is a coupling 
run, one substructure has not 
been reduced and is defined in 
BUIK. This run is a combined 
Phase 1 & 2 run. 


( REPT LOOP 37,8 ) 

* 

^Return to top of Loop 
for next substructure. 


TO LABEL L37D 
next page 
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ALTER 84 DETAIU2D FLOW 
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START GENERATING EQUILIBRIUM MATRICES 


> 

f 

NOTE: Equilibrium matrices express resultants about a chosen origin 

due to unit loads at the USET degrees of freedom. They are 
generated from an input matrix pir] and the solved structural 
transformations [Gm] , [Gc0 y and [P] . Eftr expresses 

resultants due to unit rigid body loads (r-set). 

■ 

1 

r 

PURGE EQc 

j/OMIT/EQni/MPCFl 


1 RESULTANTS /UNIT -USET LOADS | 

R 

EQ{) = [Sir] [d] 
ESULTANTS/UNIT a-S 

T 

SET LOADS 



RESULTANTS/UNIT o-SET LOADS 

. NO / OMIT ^ 

[Edo] .[Bte] [Co] 1 
RESULTANTS AJN1T f-SET LOADS 

V 7 , 

1 YES (Ni 

jj%) <= [.EQo ! EQal 

f 0 “\ 

EQUIV - - 

(kfQ* Ika] 

merged using <0?^ > 

* — 


1 > 


1 RESULTANTS /UNIT n-SET 

LOADS 


[«3 

| <= [ os ! mf) 

merged using ^ 

bih?} 


Note: Resultants due to SPC unit loads cannot 

be obtained by this method. Therefore reserve 
SPC 1 s for zero stiffness D.O.F, and sym. or 
anticon3traints at the plane of symetry. If 
there is a plane of symmetry, the resultants 
expressed should be sym. or antiresultants only. 


LABEL L8UA 



YES (No s-set) 


(fqtT 

— 

EQUIV 





f TO LABEL L81*B ’ 
V next page 
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resultants/unit 

m-SET LOADS 

= [eq/} 

[>] T 

RESULTANTS /UNIT 

g-SET LOADS 

$ 1 

ft 

f 

! EQi0 

merged using 

M 


LABEL L84B 


'MFCF1; 

Vo/ 


YES (No m-set) 


LABEL L 84 C 

\r 

MOOT SW-, 03a, BQo, BJf, B3m, B3n, BQgj 


MATCPR 


Print to check input matrix Eftr . Printed out 
in transposed form for ease of reading. 


RIGID BODY STIFFNESS 
[ xl = [krrl + furl 


MATGPR 

rxi 


— Recalculated because it wasn't 
saved in RBMG3* In case rigid 
body check doesn't worjj out, 
ve want to print out £x j . 


RESULTANTS /RIGID BOOT A'S 

D°Q -WW 

should be zero 


If not zero and If caused by an 
Isolated error, the resultants 
can help pinpoint the trouble spot. 


1 MATQPP I 

L 

LSJ 

U 


next, pace 











, L_ 

( EHD OF ALTER 8U ) 
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,• <rg* 


pppaTTKD FLOW 


ALTER 89,89 


Restrained Free modes will be 
obtained, which will probably 
be fixed at Interface (c-set). 
A col. partition matrix {CPaJ} 
must be input* 


/is\ 

RHODE 


PURGE 

Matrices generated in 
restrained mode path 


LABEL L89 


Free-Free Modes 
will be obtained 


PMfPlTIOIi REDUCED STIFF & MASS INTO 
| CONSTRAINED & FREE SET USDK5 COL. 
PARTITION MATRIX CPAJC 


CHANGE MI TO KM IN ORIGINAL 
READ STATEMENT MI HAS BEEN 
USED BIgEWHEBB 


OF ALTER 


OBTAIN RESTRAINED FREE HOPES 


HEAD KJJ, MJJ,„ »„ CASECC/LAMJ, ^/ C » * 

MODES/V, N, KEIGV ■— — 


SAVE NEIGV 

KftiTV LAMJ, LAMa/THUE 

SSlT UK W, MJ , OKOSJ 

OFP LAMJ, OEIOSJ„„/A* b » CABI3B0 


M ; 1 1 


EXPAHD EKUaiVBCTOBS "-TO A-LIKBUP 
USWJ CPAJC r -I 


JDW C*WV - -1 

0«i^<- f-™ 1 ] 

BQjUIV 

(rain)* — raija 


CHKPNT 


PHIJa 

PHIa 

LAMa 


JUMP TO 
LABEL L93 













ALTER 112 DETAILED FLOW 


(ALTER 112^ - 


This is a coupling 
run. 



YES 


PURGE 

Matrices generated 
for coupling run 


NO 


I REWIND TAPE CONTAINING REDUCED SUBSTRUCTURE MATRICES 
{ ~ If, -1/ C, N, 9/V, Y, TPNUffi9 


INITIAL LOO? VP 

JIAMETER VALUE 

wFl 

: 1 


Top of Loop to extract out -\+ JUMP L0GP112 
eigenvectors for each LABEL L00P112 


I 


( F INIS ) 



jum/ to \ 

LABEL Lila t 
next page J 


' TO LABEL L112A 
y next page 


i 

| 

i 
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REPRODUriBH ITY OF THE ORIGINAL PAGE IS POOR. 


N A S T R A N F X E C U T l V E CONTROL DECK E C H* □ 



lO PHASE? 0° l iT S'? \ 

APP D I S=> 

TIME 3? 

CHKONT T FS 

SUL 3 * 0 

PIAG 7. 13* 14. ! S* ?\ • *>? 

ALTER 2 * "* A ? AM E TER DF ’AULTS 

PAP AM / /Ct NtNO^/V . YtNOSlJM = l 

PAP AM //CfNtNP/ViY. TPCri^Y=- 1 
PAP AM //Cf N* N1P/V • Y • TPrciPYN^- 1 

PAP AM //C.N,NlP/V*Y.>? w nDF=-l 

PAP AM / /C . N t N 1P/V*N*TP UF* =— 1 

f 

ALTEP 2J f ?l * DETERMINES COUPLING PARAMETER f, PERFORMS PURGES 

PAP AM //C»N» SUS/V.N.COUOLF/Vt Y . NO SUP /C iN»? ?» 

PARAM //CtN. A>lO/VtNtNPF( mj/ V, N , NOGF NL / V • N . NCSI MP * 

PAP AM //CtNtAND/V*N,\iniNP/VtNt NOE L M T/ Y • N.CCUPL E % 

COMO EH?3tN7lN? X WILL EXIT NOT A COUPLING PUN ^ NO ELEMENTS 

PURGE G?iT 9 K GGX/N )SH° 

ChKPNT G^ST.KGGX 

CPNP L ? 3 1 COUPLE 

PURGE OG^ST/TIU* 

JUMP L2A 

LAHTL ER?3 

PWTPAPM //C.N» o/c.NtNniNP & 

FX IT 

LAHFL L2J 

Y 

At TEP ?A 
LABEL L27 

At TER 25 

COND L2 7tNOSlM? 

AITFR 2 7 
LABEL L 27 

ALTER ?9 

C CNO L >Rf COUPLE 

JUMP L 3 0 

L Af' l ’l L 29 

ALTER 70 
LABEL L 30 

I 

AITFR 77 * * 7 i |F COUPLING RUM .COMBI NFS SUF» S TPUC TURES • 

* THE PILLOWING MATPICFS FOP EACH SUOSTPUCT. ASSUMED ON IN 

$ C? G I - COLUMN PARTI TIONG VECTOR F OR MERGING Kl & Ml 

A < I r * Ml = REDUCED STI^FKFISS & MASS 

PURGE C?GI .< It M U<GGT .MGS I .KGGS. MSG S t KG T t MG T/C CU PLE 

COND L ** ?O f C0li°L ~ % SKI" .NOT A COUPLING RUN 

U’PUTTl /» t « t^CtNt - 3/C .NtS/Vt Yf T fl NA Mr g * LIST TA PF f. REWIND 

PARAM //CtN.N r )*/VtN t n ASS*t Y INITIAL LOOP PASS PARAMETER 
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K 





** ' • 


‘ ; REPRODUCIBII ITY OF THE ORIGINAL PAGE IS POOR. 


N A S T R A N| 


F X p C U T I V F 


C 0 fs T R C L 


DECK 


ECHO 


JUMP Lm 3 ^7 

I ABEL LOOP 37 S T )P OF 

[PAP AM / /C f M* SUB/V. N*PA 3S 1 / V »N »PA S S/C , H , 2 

| INPUTT 1 /C’GI.Kl.MT, , /C.N.O/C.N.Q * 

jCONO L37 f ^ASSl * SK 1 3 TO L37 IF FIRST PASS 

[jump L3 7A 

juARFL L V 

[MERGE. . . • < I . C 3 G I • /KGGS/C . N . - 1 /C . N , 2 /C . N . ft 

[MERGE. • . , M I.C^GI • /MGGS/C tN.-l/C.N.P/C.N.f* 

^L ABEL L 3 7 A 

COND L S SKIP TO L37B TF FIRST PASS 

MFRGF, . . .< I. C >GI . /KGr, | /C ,N 1 /C .N .2 /C .N.ft 

m FP GF 4 . . . M I . C 3 j t • /MGG t /C . N . - 1 /C .N*?/C .N.6 

ADD < GG3.<GSl/*GT S 

FOJ TV K GT. <GGS/TP1JF 

ADD MGGS.MGGI /MGT S 

FOJ IV M GT. MGGS'TRJE 

LAUFL L 37 I 

PAR AM //O.N. AO) /V.N.R ASS/ V.\l. PASS/C .N.l 

PAR AM / /C. N. SUO/V. M. SK t P2/V * Y . NO SUB/V * N. PASS 

CONP L17C.NJFL4T 

PARAM //CtN*StJT/V.NtSKIR?/V.NfSKTP2/CfN.l 

LAPEL L }7Z 

CP NO L J7 O, S< IR > 

PFPT L r) jp 37, \ 

LABEL L 370 

fHKPNT < GGS.MGGS 

A 00 < GGX ,K GGS/<GGY * 

AOD MCG. MGG3/4GGY * 

CHKPNT < G GY . M GGY 

FOUIV < GGV ♦ < GG/ vIOGTNL T 


altep 
sm ab 

At. T f R 
ADO 

MAT GPP 

FOU IV 

ALTtP 

CONO 

JUMP 

LABEL 

alter 

MCE 2 
AL TER 
SEEM AT 
$ 

A! TfcP 
PURGE 
PURGE 


40, *0 

CE I. <GG7/<GG/V,N,LUS Jr T/V ,M ,NOGE NL/ V . N . NOS I MP 
48, AS 

MGS. /VGG/C./,ALPHA=( T«^.4,0. O) % 

GPL . U S F T t sit t V/GG//C * N * G 
<GG.KNN/MPC c l/MGGY.MNN/MPCFl 
40 

LA >. COUPLE 
L 01.4 
LAO 
E*. *3® 

U SET , GM . K GG . MGG Y . ,/KNK.MNN*, * 

74 

K A M A A « . . //C * N . PR I M T 


CP V<L. OA . CPNSF . rPGMM .FUR .E OL .E 0 A . E OH . E QF ,E ON • E QM , FQG/ R E ACT 

fX . EXT, f )MT. ^ONT.EOGT.f 3GTC * MHGCi . MOGG Y/Rf ACT 
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N A 3 T D % M 


r u t i v f 


c P k t « r: L 


0 F C K 


F C 


CONO 

PViPGF 

VFC 

CONO 

V FT 

LA^FL 

COND 

VFC 

I APFL 

CONO 

VFC 

LAHR 

CH< PNT 

IUJPGF 

TRNSP 

MPY AD 

MF9GF 

FOJ IV 

CONO 

tonsp 
MPY AD 
MFPGF 
L AOFL 
FOIJ IV 
CONO 
MFPGF, 
LARFL 
COU IV 
COND 
TPNSP 
MPY AD 

mepge 


L94F*OFACT T MtGT 1 F l)P r INFO Tn GF NF u A T F FQG 

C^FOA/OM I F/C^NSF/Sl MGL P /CPGMN/MPC- 1 
U 3FT /C 3 A / C • N • A /C t N * F /C , N * L * 

LC°ltlM IT 

u < iFT/C > ' r ^/ r . «m »F/r *N-«n/r *n*a % 

L C» t 

L 0^2 •SINGLE 

lis FT /C*» 4SF/C .N- t N/C |N* C ./C»N»P * 

L C* 3 ? 

L C r> 3 t M >CF I 

1 1 S E T /C GM 4 /0 * N * G /C * N * M /C * N « N ft 

L C r > 3 
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TABLE 1 .- STATISTICAL DESCRIPTION OF l/8-SCALE MODEL 
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Schematic of cargo door joints. 





















A SUBSTRUCTURES 
MUST HAY*! UNIQUE 
GRID I.D. NUMBERS 


OUTPUT OK TARE 

REDUCED PSEUDO -STRUCTURE 
STI FFN E SS A MASS Ka^ # Maa 
FOR POSSIBLE FURTHER 
COUPLING. 


A TAPE GENERATED 
FOR EACH SUBSTRUCTURE 
CONTAINING RESULTS FROM 
PHASE 2 SO THAT DETAIIED 
INFORMATION CAN BE OBTAINED 
IN PHASE 3 


NASTRAN 
PHASE 3 


Figure 21. • Flow diagram for NASTRAN subetructuring to obtain normal 

modes ( Rigid Format J). 
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STA. 118.16 


UNDEFORMED 

DEFORMED 

( FREQUENCY ■ 56 . U Hz) 



(a) First and second free bending modes. 
Figure 30.- Shape for SRB modes. 







INNER RADIUS OP 
PROPELLANT 



CENTER LINE 
OF SRB 


(b) Free longitudinal rod mode showing longitudinal thickness 
shear deflection. I 96 .O Hz. 


Figure 30.- Concluded. 
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(a) First free bending mode. 
56.4 Hz. 


(b) Second free bending (c) Longitudinal mode 
mode. 173.0 Hz. showing some tor- 

sion. 196.1 Hz. 
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Figure 31.- Shapes for SRB bending modes. 


